Forecasting effects of climate changes on the population genetic structure of Anacardium occidentale in the Cerrado biome, Brazil.

There has been a continuous interest in understanding the patterns of genetic diversity in natural populations because of the role of intraspecific genetic diversity as the basis of all evolutionary change and thus, its potential effects on population persistence when facing environmental changes. Here, we provided the first description of genetic diversity distribution and population structure of Anacardium occidentale L. (cashew) from the Brazilian Cerrado, one of the most economically important tropical crops in the world. We applied Bayesian clustering approaches (STRUCTURE and POPS) that allow predicting the effects of future climatic changes on the population genetic structure of A. occidentale. We identified distinct genetic groups corresponding to the southwestern, central, and northern regions of the species' range. The characterized genetic clusters will disappear under future climate change scenarios, leading to a homogenization of genetic variability across the landscape. Our findings suggest a high likelihood for the loss of genetic diversity, which in turn will reduce the evolutionary potential of the species to cope with predicted future climatic changes. Results from this study may help develop management strategies to conserve the genetic diversity and structure of natural cashew populations.

Holocene climate changes explain the spatial pattern in genetic diversity in populations of Cyperus papyrus from Southeast Africa wetlands.

Wetlands are one of the most threatened ecosystems in the world because more than 70% of the area worldwide has been lost since 1900. Wetland plant species rely greatly on water for seeds and propagules, which may lead to a downstream unidirectional dispersal and accumulation of genetic diversity downstream. However, several species show no support for unidirectional genetic diversity, revealing the complexity of population dynamics and gene flow in wetlands. Here, we used microsatellite loci to address how the past demographic dynamics shaped the contemporary spatial pattern in genetic diversity and population structure of Cyperus papyrus in wetlands of Southeast Africa. Using spatially explicit analysis and coalescent modelling, we found no support for unidirectional dispersal. Instead, we found higher genetic diversity in populations upstream than downstream in the river basin. We also found high admixture among populations, most likely due to connections between adjacent river basins during sporadic floods, and ongoing gene flow due to bird-mediated seed dispersal. Our results suggest stepping-stone migration due to strong isolation-by-distance, but not necessarily unidirectional. Moreover, the past demographic dynamics in the Holocene shaped the current pattern of genetic diversity and structure, leading to higher genetic diversity in populations upstream the Zambezi river basin. Our results also point to the very low genetic diversity of C. papyrus populations in Southeast Africa and the need for management and conservation strategies to guarantee the long-term persistence of the species in the region.

Chloroplast genome assembly of Handroanthus impetiginosus: comparative analysis and molecular evolution in Bignoniaceae.

MAIN CONCLUSION: Bignoniaceae species have conserved chloroplast structure, with hotspots of nucleotide diversity. Several genes are under positive selection, and can be targets for evolutionary studies. Bignoniaceae is one of the most species-rich family of woody plants in Neotropical seasonally dry forests. Here we report the assembly of Handroanthus impetiginosus chloroplast genome and evolutionary comparative analyses of ten Bignoniaceae species representing the genera for which whole-genome chloroplast sequences were available. The chloroplast genome of H. impetiginosus is 159,462 bp in size and has a similar structure compared to the other nine species. The total number of genes was slightly variable amongst the Bignoniaceae, ranging from 124 in H. impetiginosus to 144 in Anemopaegma acutifolium. The inverted repeat (IR) size was variable, ranging from 24,657 bp (Tecomaria capensis) to 40,481 bp (A. acutifolium), due to the contraction and retraction at its boundaries. However, gene boundaries were very similar among the ten species. We found 98 forward and palindromic dispersed repeats, and 85 simple sequence repeats (SSRs). In general, chloroplast sequences were highly conserved, with few nucleotide diversity hotspots in the genes accD, clpP, rpoA, ycf1, ycf2. The phylogenetic analysis based on 77 coding genes was highly consistent with Angiosperm Phylogeny Group (APG) IV. Our results also indicate that most genes are under negative selection or neutral evolution. We found no evidence of branch-site selection, implying that H. impetiginosus is not evolving faster than the other species analyzed, notwithstanding we found site positive selection signal in several genes. These genes can provide targets for evolutionary studies in Bignoniaceae and Lamiales species.

Multi-Scale Landscape Influences on Genetic Diversity and Adaptive Traits in a Neotropical Savanna Tree.

Changes in landscape structure can affect essential population ecological features, such as dispersal and recruitment, and thus genetic processes. Here, we analyze the effects of landscape metrics on adaptive quantitative traits variation, evolutionary potential, and on neutral genetic diversity in populations of the Neotropical savanna tree Tabebuia aurea. Using a multi-scale approach, we sampled five landscapes with two sites of savanna in each. To obtain neutral genetic variation, we genotyped 60 adult individuals from each site using 10 microsatellite loci. We measured seed size and mass. Seeds were grown in nursery in completely randomized experimental design and 17 traits were measured in seedlings to obtain the average, additive genetic variance (V a ) and coefficient of variation (CV a %), which measures evolvability, for each trait. We found that habitat loss increased genetic diversity (He) and allelic richness (AR), and decreased genetic differentiation among populations (F ST ), most likely due to longer dispersal distance of pollen in landscapes with lower density of flowering individuals. Habitat amount positively influenced seed size. Seeds of T. aurea are wind-dispersed and larger seeds may be dispersed to short distance, increasing genetic differentiation and decreasing genetic diversity and allelic richness. Evolvability (CV a %) in root length decreased with habitat amount. Savanna trees have higher root than shoot growth rate in the initial stages, allowing seedlings to obtain water from water tables. Landscapes with lower habitat amount may be more stressful for plant species, due to the lower plant density, edge effects and the negative impacts of agroecosystems. In these landscapes, larger roots may provide higher ability to obtain water, increasing survival and avoiding dying back because of fire. Despite the very recent agriculture expansion in Central Brazil, landscape changes are affecting neutral and adaptive variation in T. aurea. Several populations have low additive genetic variation for some traits and thus, may have limited evolvability, which may jeopardize species long-term persistence. The effect of habitat loss on highly variable neutral loci may only be detected after a certain threshold of population size is attained, that could become dangerously small masking important losses of heterozygosity endangering species conservation.

Agricultural Landscape Heterogeneity Matter: Responses of Neutral Genetic Diversity and Adaptive Traits in a Neotropical Savanna Tree.

Plants are one of the most vulnerable groups to fragmentation and habitat loss, that may affect community richness, abundance, functional traits, and genetic diversity. Here, we address the effects of landscape features on adaptive quantitative traits and evolutionary potential, and on neutral genetic diversity in populations of the Neotropical savanna tree Caryocar brasiliense. We sampled adults and juveniles in 10 savanna remnants within five landscapes. To obtain neutral genetic variation, we genotyped all individuals from each site using nine microsatellite loci. For adaptive traits we measured seed size and mass and grown seeds in nursery in completely randomized experimental design. We obtained mean, additive genetic variance (V a ) and coefficient of variation (CV a %), which measures evolvability, for 17 traits in seedlings. We found that landscapes with higher compositional heterogeneity (SHDI) had lower evolutionary potential (CV a %) in leaf length (LL) and lower aboveground dry mass (ADM) genetic differentiation (Q ST ). We also found that landscapes with higher SHDI had higher genetic diversity (He) and allelic richness (AR) in adults, and lower genetic differentiation (F ST ). In juveniles, SHDI was also positively related to AR. These results are most likely due to longer dispersal distance of pollen in landscapes with lower density of flowering individuals. Agricultural landscapes with low quality mosaic may be more stressful for plant species, due to the lower habitat cover (%), higher cover of monocropping (%) and other land covers, and edge effects. However, in landscapes with higher SHDI with high quality mosaic, forest nearby savanna habitat and the other environments may facilitate the movement or provide additional habitat and resources for seed disperses and pollinators, increasing gene flow and genetic diversity. Finally, despite the very recent agriculture expansion in Central Brazil, we found no time lag in response to habitat loss, because both adults and juveniles were affected by landscape changes.

Diversification of the widespread neotropical frog Physalaemus cuvieri in response to Neogene-Quaternary geological events and climate dynamics.

Here we reconstructed the demographical history and the dispersal dynamics of Physalemus cuvieri through the Neogene-Quaternary periods by coupling DNA regions with different mutation rates, ecological niche modelling, reconstruction of spatio-temporal lineage dispersal and coalescent simulations. Still, to test alternative diversification scenarios we used approximate Bayesian computation. Molecular phylogenetic analysis recovered four deep and strongly supported clades, which we interpret as population lineages. The ancestral location reconstruction placed the root in southcentral Amazonia, and the dispersal events indicate that spatial displacement was widespread early in the diversification of this species. The demographical scenario of "Multiple Refugia" with recent lineage admixture was the most likely hypothesis to predict the observed genetic parameters of P. cuvieri. Our results revealed that Neogene orogenic events might have played a prominent role in the early diversification of P. cuvieri. The species shows deep divergences with strong regional population structure, despite its widespread distribution. Final uplift of the central Brazilian Plateau and formation of the river basins in Central South America played an important role in the origin, diversification and the maintenance of P. cuvieri lineages.

First description and evaluation of SNPs in the ADH and ALDH genes in a population of alcoholics in Central-West Brazil.

Worldwide, different studies have reported an association of alcohol-use disorder (AUD) with different types of Single Nucleotide Polymorphisms (SNPs) in the genes for aldehyde dehydrogenase (ALDH) and alcohol dehydrogenase (ADH). In Brazil, there is little information about the occurrence of these SNPs in the AUD population and an absence of studies characterizing the population in the Central-West Region of Brazil. Actually, in Brazil, there are more than 4 million people with AUD. Despite the major health hazards of AUD, information on alcohol consumption and its consequences are not well understood. Therefore, it is extremely important to characterize these SNPs for the better understanding of AUD as a genetic disease in the Brazilian population. The present study, unlike other studies in other countries, is done with a subject population that shows a significant amount of racial homogenization. We evaluated the presence of SNPs in the ADH (ADH1B, ADH1C, and ADH4) and ALDH (ALDH2) genes in alcohol users of Goiânia, State of Goiás - Brazil, and then we established a possible relationship with AUD by allelic and genotypic study. This study was conducted with a population of people with AUD (n = 99) from Goiás Alcohol Dependence Recovery Center (GO CEREA) and Psychosocial Care Center for Alcohol and Drugs (CAPS AD), and with a population of people without AUD as controls (n = 100). DNA was extracted from whole-blood samples and the genotyping was performed using TaqMan® SNP genotyping assays. For characterization and evaluation of SNPs in the population, genotype frequency, allele frequency, haplotype frequency, Hardy-Weinberg equilibrium, and linkage disequilibrium were analyzed. Statistical analyses were calculated by GENEPOP 4.5 and Haploview software. The allele 1 was considered as "wild" (or *1) and allele 2 as mutant (or *2). Significant differences were found for ADH1B*, ADH4*2, and ALDH2*2 SNPs when the genotype and allele frequencies were analyzed. In addition, four haplotypes were observed between ADH1B*2 and ADH1C*2 through linkage disequilibrium analysis. The genetic variants may be associated with protection against AUD in the population studied.

A large historical refugium explains spatial patterns of genetic diversity in a Neotropical savanna tree species.

BACKGROUND AND AIMS: The relative role of Pleistocene climate changes in driving the geographic distribution and genetic diversity of South American species is not well known, especially from open biomes such as the Cerrado, the most diverse tropical savanna, encompassing high levels of endemism. Here the effects of Quaternary climatic changes on demographic history, distribution dynamics and genetic diversity of Dimorphandra mollis, an endemic tree species widely distributed in the Cerrado, were investigated. METHODS: A total of 38 populations covering most of the distribution of D. mollis were analysed using internal transcribed spacer (ITS) sequences and nuclear microsatellite variation [simple sequence repeats (SSRs)]. The framework incorporated statistical phylogeography, coalescent analyses and ecological niche modelling (ENM). KEY RESULTS: Different signatures of Quaternary climatic changes were found for ITS sequences and SSRs corresponding to different time slices. Coalescent analyses revealed large and constant effective population sizes, with high historical connectivity among the populations for ITS sequences and low effective population sizes and gene flow with recent population retraction for SSRs. ENMs indicated a slight geographical range retraction during the Last Glacial Maximum. A large historical refugium across central Brazil was predicted. Spatially explicit analyses showed a spatial cline pattern in genetic diversity related to the paleodistribution of D. mollis and to the centre of its historical refugium. CONCLUSIONS: The complex genetic patterns found in D. mollis are the result of a slight geographical range retraction during the Last Glacial Maximum followed by population expansion to the east and south from a large refugium in the central part of the Cerrado. This historical refugium is coincident with an area predicted to be climatically stable under future climate scenarios. The identified refugium should be given high priority in conservation polices to safeguard the evolutionary potential of the species under predicted future climatic changes.

Defence responses in rice plants in prior and simultaneous applications of Cladosporium sp. during leaf blast suppression.

An alternative method to control rice blast (Magnaporthe oryzae) is to include biological agent in the disease management strategy. The objective of this study was to assess the leaf blast-suppressing effects of rice phylloplane fungi. One Cladosporium sp. phylloplane fungus was shown to possess biocontrolling traits based on its morphological characteristics and an analysis of its 18S ribosomal DNA. Experiments aimed at determining the optimal time to apply the bioagent and the mechanisms involved in its rice blast-suppressing activities were performed under controlled greenhouse conditions. We used foliar spraying to apply the Cladosporium sp. 48 h prior to applying the pathogen, and we found that this increased the enzymatic activity. Furthermore, in vitro tests performed using isolate C24 showed that it possessed the ability to secrete endoxylanases and endoglucanases. When Cladosporium sp. was applied either prior to or simultaneous with the pathogen, we observed a significant increase in defence enzyme activity, and rice blast was suppressed by 84.0 and 78.6 %, respectively. However, some enzymes showed higher activity at 24 h while others did so at 48 h after the challenge inoculation. Cladosporium sp. is a biological agent that is capable of suppressing rice leaf blast by activating biochemical defence mechanisms in rice plants. It is highly adapted to natural field conditions and should be included in further studies aimed at developing strategies to support ecologically sustainable disease management and reduce environmental pollution by the judicious use of fungicidal sprays.

Ancestral reconstruction of reproductive traits shows no tendency toward terrestriality in leptodactyline frogs.

BACKGROUND: Traditionally, the evolution of terrestrial reproduction in anurans from ancestors that bred in water has been accepted in the literature. Still, the existence of intermediate stages of water dependency, such as species that lay eggs close to water (e.g., in burrows) instead of in bodies of water, supports the hypothesis of an ordered and gradual evolution in the direction of a more terrestrial form of reproduction. However, this conventional view has recently been challenged for some anurans groups. Leptodactylinae frogs are a remarkable example of anurans with an outstanding diversity in terms of reproductive features, with distinct water dependency among lineages. Here, we tested the hypothesis of a gradual and ordered tendency towards terrestriality in Leptodactylinae, including the existence of obligatory intermediate stages, such as semi-terrestrial reproductive strategies. We also addressed the association between reproductive modes and the morphological and ecological features. RESULTS: An ancestral reconstruction analysis indicated that even though shifts from aquatic to terrestrial breeding occurred throughout the history of Leptodactylus and Adenomera, shifts from terrestrial to aquatic reproduction happened at almost the same frequency. Our results also demonstrated that reproductive modes for semi-terrestrial tadpoles were not necessarily an intermediate form between aquatic and terrestrial breeds. Correlations among reproductive modes and other life-history traits suggested that tadpole environment, clutch size, nuptial spines, and egg pigmentation were co-evolving and driven by water dependency. CONCLUSIONS: Our results found no evidence of evolutionary tendencies toward terrestriality in Leptodactylinae. We found reversals from terrestrial to aquatic tadpole development and no evidence of obligatory intermediate stages, such as semi-terrestrial reproductive strategies. We also found correlations between reproductive modes and other life-history traits driven by water dependence. Aquatic reproductive modes are associated with higher clutch sizes, lentic waters, and the presence of nuptial spines and egg pigmentation.

A new eigenfunction spatial analysis describing population genetic structure.

Several methods of spatial analyses have been proposed to infer the relative importance of evolutionary processes on genetic population structure. Here we show how a new eigenfunction spatial analysis can be used to model spatial patterns in genetic data. Considering a sample of n local populations, the method starts by modeling the response variable (allele frequencies or phenotypic variation) against the eigenvectors sequentially extracted from a geographic distance matrix (n × n). The relationship between the coefficient of determination (R(2)) of the models and the cumulative eigenvalues, which we named the spatial signal-representation (SSR) curve, can be more efficient than Moran's I correlograms in describing different patterns. The SSR curve was also applied to simulated data (under distinct scenarios of population differentiation) and to analyze spatial patterns in alleles from microsatellite data for 25 local populations of Dipteryx alata, a tree species endemic to the Brazilian Cerrado. The SSR curves are consistent with previous phylogeographical patterns of the species, revealing combined effects of isolation-by-distance and range expansion. Our analyses demonstrate that the SSR curve is a useful exploratory tool for describing spatial patterns of genetic variability and for selecting spatial eigenvectors for models aiming to explain spatial responses to environmental variables and landscape features.

Geographical patterns of turnover and nestedness-resultant components of allelic diversity among populations.

The analysis of geographical patterns in population divergence has always been a powerful way to infer microevolutionary processes involved in population differentiation, and several approaches have been used to investigate such patterns. Most frequently, multivariate spatial patterns of population differentiation are analyzed by computing pairwise genetic distances or F(ST) (or related statistics, such as ϕ(ST) from AMOVA), which are then correlated with geographical distances or landscape features. However, when calculating distances, especially based on presence-absence of alleles in local populations, there would be a confounding effect of allelic richness differences in the population differentiation. Moreover, the relative magnitude of these components and their spatial patterns can help identifying microevolutionary processes driving population differentiation. Here we show how recent methodological advances in ecological community analyses that allows partitioning dissimilarity into turnover (turnover) and richness differences, or nestedness-resultant dissimilarity, can be applied to allelic variation data, using an endemic Cerrado tree (Dipteryx alata) as a case study. Individuals from 15 local populations were genotyped for eight microsatellite loci, and pairwise dissimilarities were computed based on presence-absence of alleles. The turnover of alleles among populations represented 69 % of variation in dissimilarity, but only the richness difference component shows a clear spatial structure, appearing as a westward decrease of allelic richness. We show that decoupling richness difference and turnover components of allelic variation reveals more clearly how similarity among populations reflects geographical patterns in allelic diversity that can be interpreted in respect to historical range expansion in the species.

Phylogeography of Tibouchina papyrus (Pohl) Toledo (Melastomataceae), an endangered tree species from rocky savannas, suggests bidirectional expansion due to climate cooling in the Pleistocene.

Many endemic species present disjunct geographical distribution; therefore, they are suitable models to test hypotheses about the ecological and evolutionary mechanisms involved in the origin of disjunct distributions in these habitats. We studied the genetic structure and phylogeography of Tibouchina papyrus (Melastomataceae), endemic to rocky savannas in Central Brazil, to test hypothesis of vicariance and dispersal in the origin of the disjunct geographical distribution. We sampled 474 individuals from the three localities where the species is reported: Serra dos Pirineus, Serra Dourada, and Serra de Natividade. Analyses were based on the polymorphisms at cpDNA and on nuclear microsatellite loci. To test for vicariance and dispersal we constructed a median-joining network and performed an analysis of molecular variance (AMOVA). We also tested population bottleneck and estimated demographic parameters and time to most recent common ancestor (TMRCA) using coalescent analyses. A remarkable differentiation among populations was found. No significant effect of population expansion was detected and coalescent analyses showed a negligible gene flow among populations and an ancient coalescence time for chloroplast genome. Our results support that the disjunct distribution of T. papyrus may represent a climatic relict. With an estimated TMRCA dated from ∼836.491 ± 107.515 kyr BP (before present), we hypothesized that the disjunct distribution may be the outcome of bidirectional expansion of the geographical distribution favored by the drier and colder conditions that prevailed in much of Brazil during the Pre-Illinoian glaciation, followed by the retraction as the climate became warmer and moister.

Development of microsatellite markers for the neotropical tree species Dipteryx alata (Fabaceae).

PREMISE OF THE STUDY: Microsatellite markers were developed for the population genetic analyses of the neotropical tree Dipteryx alata (Fabaceae). METHODS AND RESULTS: Microsatellites were developed from a genomic shotgun library. Polymorphism at each microsatellite loci was analyzed based on 94 individuals from three populations. Eight loci amplified successfully and presented one to 10 alleles, and expected heterozygosities ranged from 0.097 to 0.862. Four loci also amplified in Pterodon emarginatus and presented similar polymorphism. CONCLUSION: The eight microsatellite primer pairs are potentially suitable for population genetic studies and successfully amplified in another Fabaceae species.

A geographical genetics framework for inferring homing reproductive behavior in fishes.

One of the most intriguing patterns of migration and gene flow that affects genetic structure is the reproductive homing behavior of fishes, wherein the adults return to the areas in which they were spawned. Here we reviewed the literature on homing behavior in fish and propose an analytical framework for testing hypotheses regarding this behavior and its effects on the genetic structure of fish in an explicit geographical context, using a geographical genetics toolbox. Although disentangling the many potential causes underlying genetic population structure and unambiguously demonstrating that the homing behavior causes these genetic patterns is difficult, our framework allows the successive testing of homing behavior with increasing levels of complexity based on the following: (1) establishment of population structures among waterheads; (2) patterns of genetic variability throughout the adult migratory pool; (3) analyses of the non-migratory adult pool; and (4) comparisons among successive generations. We expect that the framework presented here will help delineating the appropriate uses of different sampling designs to make inferences regarding homing behavior and illustrate the limits imposed by the interpretation of different types of genetic data. More importantly, we hope this framework enables researchers to understand how a particular dataset can be utilized in a broader context as an ongoing part of a larger research program and thus guide future research by developing better and more integrated sampling designs.

Late Miocene diversification and phylogenetic relationships of the huge toads in the Rhinella marina (Linnaeus, 1758) species group (Anura: Bufonidae).

We investigated the phylogeny and biogeography of the Rhinella marina group, using molecular, morphological, and skin-secretion data, contributing to an understanding of Neotropical faunal diversification. The maximum-parsimony and Bayesian analyzes of the combined data recovered a monophyletic R. marina group. Molecular dating based on Bayesian inferences and fossil calibration placed the earliest phylogenetic split within the R. marina group at ∼ 10.47 MYA, in the late Miocene. Two rapid major diversifications occurred from Central Brazil, first northward (∼ 8.08 MYA) in late Miocene and later southward (∼ 5.17 MYA) in early Pliocene. These results suggest that barriers and dispersal routes created by the uplift of Brazilian Central Shield and climatic changes explain the diversification and current species distributions of the R. marina group. Dispersal-vicariance analyzes (DIVA) indicated that the two major diversifications of the R. marina group were due to vicariance, although eleven dispersals subsequently occurred.

Phylogeography and disjunct distribution in Lychnophora ericoides (Asteraceae), an endangered cerrado shrub species.

BACKGROUND AND AIMS: Lychnophora ericoides (Asteraceae) presents disjunct geographical distribution in cerrado rupestre in the south-east and central Brazil. The phylogeography of the species was investigated to understand the origin of the disjunct geographical distribution. METHODS: Populations in the south and centre of Serra do Espinhaço, south-east Brazil and on ten other localities in Federal District and Goiás in central Brazil were sampled. Analyses were based on the polymorphisms at chloroplast (trnL intron and psbA-trnH intergenic spacer) and nuclear (ITS nrDNA) genomes. From 12 populations, 192 individuals were sequenced. Network analysis, AMOVA and the Mantel test were performed to understand the relationships among haplotypes and population genetic structure. To understand better the origin of disjunct distribution, demographic parameters and time to most recent common ancestor (T(MRCA)) were estimated using coalescent analyses. KEY RESULTS: A remarkable differentiation between populations from the south-east and central Brazil was found and no haplotype was shared between these two regions. No significant effect of isolation by distance was detected. Coalescent analyses showed that some populations are shrinking and others are expanding and that gene flow between populations from the south-east and central Brazil was probably negligible. CONCLUSIONS: The results strongly support that the disjunct distribution of L. ericoides may represent a climatic relict and that long-distance gene flow is unlikely. With an estimated time to most recent common ancestor (T(MRCA)) dated from approx. 790,655 +/- 36,551 years bp (chloroplast) and approx. 623,555 +/- 55,769 years bp (ITS), it was hypothesized that the disjunct distribution may be a consequence of an expansion of the geographical distribution favoured by the drier and colder conditions that prevailed in much of Brazil during the Kansan glaciation, followed by the retraction of the distribution due to the extinction of populations in some areas as climate became warmer and moister.

Molecular dating and phylogenetic relationships among Teiidae (Squamata) inferred by molecular and morphological data.

We present a phylogenetic analysis of teiid lizards based on partitioned and combined analyses of 12S and 16S mitochondrial DNA sequences, and morphological and ultrastructural characters. There were some divergences between 12S and 16S cladograms, but phylogenetic analyses of the combined molecular data corroborated the monophyly of Tupinambinae, Teiinae, and "cnemidophorines", with high support values. The total combined analysis (molecules+morphology) produced similar results, with well-supported Teiinae and "cnemidophorines". We present an evolutionary scenario for the evolution of Teiidae, based on molecular dating of evolutionary events using Bayesian methods, ancestral areas analysis, the fossil record, the geographic distribution of genera, and environmental and geologic changes during the Tertiary. According to this scenario, (1) all current teiid genera, except Aspidoscelis, originated in isolation in South America; (2) most teiid genera originated during the Eocene, a period characterized by savanna expansion in South America; and (3) Cnemidophorus originated in South America, after which some populations dispersed to Central America during the Late Miocene.